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o 

i We calculate the induced elastic-interaction between pointwise membrane inclusions that locally 

, interact up to quadratic order with the membrane curvature tensor. For isotropic inclusions, we 

recover the usual interaction proportional to the inverse fourth power of the separation, however 
with a prefactor showing a non-trivial dependence on the rigidity P of the quadratic potential. In the 
large F limit, corresponding to "hard" inclusions, we recover the standard prefactor first obtained 
by Goulian et al. [Europhys. Lett. 22, 145 (1993)]. In the small T limit, corresponding to "soft" 
inclusions, we recover the recent result of Marchenko and Misbah [Eur. Phys. J. E 8, 477 (2002)]. 
This shows that the latter result bears no fundamental discrepancy with previous works, but simply 
corresponds to the limit of soft inclusions. We discuss how the same inclusion can be depicted as 
hard or soft according to the degree of coarse-graining of the pointwise description. Finally, we 
argue that conical transmembrane proteins should be fundamentally considered as hard inclusions. 

PACS numbers: 87.16.Dg (Membranes, bilayers, and vesicles), 87. 15. Kg (Molecular interactions; membrane- 
protein interactions) 



^ I. INTRODUCTION 

a, 

In a recent paper [J, Marchenko and Misbah have revisited the calculation of the long-range, membrane- mediated 
i^i - interaction acting between inclusions embedded in biological or surfactant membranes. Using a method inspired 
from electrostatics, they calculated an "induced interaction" , which they compared with the results of several works 
t-H , previously published on the subject @, 0, 0, 0, ■ They pointed out two major discrepancies: (i) their result is 
proportional to JJk 2 (k being the bending rigidity of the membrane) while the induced interaction calculated in 
Refs. HHSIllg is directly proportional to K, (ii) the result of Refs. involves the size of the inclusion 

' (either directly or through a cutoff in the case of pointwise descriptions) whereas the size of the inclusion does not 
! appear in the result of Ref. . 

In this contribution, we show that there are no discrepancies between the approaches of Refs. Q H, Q, |^ El an d 
that of Ref. Q, neither in the spirit of the model nor in the results. We demonstrate that the different expressions 
obtained for the induced interaction correspond to different limits of a common, more general model. We do so by 
solving a pointwise quadratic interaction model similar to the one used in Ref. [lj , exactly, instead of perturbatively, 
and without the need to introduce a high-wavector cutoff. According to the rigidity T of the quadratic interaction 
potential, we find either the standard prefactor first obtained by Goulian et al. @, M El El , in the large F limit 
1 , | corresponding to "hard" inclusions, or the result of Marchenko and Misbah pj, in the small T limit corresponding to 
"soft" inclusions. 

In the second part of the paper we discuss the role of the cutoff (introduced e.g. in Refs. 0,01) in relation with a 
coarse-graining procedure of the pointwise description. We determine how the effective quadratic interaction potential 
depends on the coarse-graining scale, and we discuss the distinction between hard and soft inclusions in the light of 
this analysis. 
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II. INDUCED INTERACTION 



Here, we calculate the induced interaction between pointwise "defects" , or inclusions. We describe the inclusions 
using the same formalism as in Ref. We use, however, a more general Hamiltonian to describe the elasticity of 
the membrane, and, in addition, we perform our calculations exactly, whereas the result of Ref. []] was obtained at 
lowest-order in a perturbative expansion. We emphasize that we shall not introduce any cutoff during our calculation. 

The generalized membrane elastic Hamiltonian we use is 



F m = - I dxdy 



-yh 2 + a (V/i) 2 + k (V 2 /i) 2 + k g (V V 2 /i) 2 + k 8 (V 4 /i) 



(1) 
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where h(x,y) describes the normal displacement of the membrane. Being quadratic in h(x,y), this Hamiltonian is 
accurate only for small membrane deformations. We have assumed that the membrane has no physical boundary 
except at infinity, which allowed us to discard all the terms that could be written as the divergence of a vector field 
(e.g., the Gaussian curvature energy 0). The parameters 7 and a describe a confining potential and a membrane 
tension, respectively. In most cases they are negligibly small, and they play no role except at very large scales. The 
parameter n is the usual bending constant and Ke, k$ are higher-order bending elasticity constants We may 
rewrite the membrane Hamiltonian in the reciprocal space as 



F m = 



d 2 q k 



+ cV + q 4 (1 + bU 2 + ^ 4 )]|Mq)l 2 



(2) 



Here L is the lateral extension of the membrane; be = {k,q/ k) 1 ' 2 , b$ — (ks/k) 1 ^ 4 are nanometric lengths of the order of 
the thickness of the membrane, and £ 7 = (K/7) 1 / 4 , £ a = (k/ct) 1 / 2 are very large lengths in the case of free membranes 
with a vanishing tension. 

To describe the coupling of the membrane with a pointwise isotropic inclusion labeled i, situated in r = r^, we 
generalize the expression used in Ref. [lj (in which the term proportional to D below was omitted): 



U, 



dx dy 6(r — 



AV 2 h{v) 



D 



d a dph(T) d a d/3h(i 



D 



d a d a h(r) dpdphiY) 



(3) 



where q,/9 £ {1, 2} stand for the x and y coordinates and 8{r) is the Dirac distribution. Summation over repeated 
indices will be implied throughout. For a membrane containing two isotropic inclusions, the total energy functional 
is F m + Ui + U2- To proceed, we rewrite the interaction terms Ui, up to a constant, as 



Ui = 



(0 

a/3 



c. 
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(4) 



where H 



(0 - 



d a dph{vi), and C^p = cq 5 a p for an isotropic inclusion. The constants T, e and cq are given by Y = —D, 



L a{3 

€ = D/D and cq = A/D. Note that the condition e > — i has to be satisfied, overwise the membrane would be locally 
unstable. In order to give a simple interpretation of the interaction, we can also write Eq. (0J in term of the principal 
values and H% of the local curvature tensor H^l, i.e., 



CT i = -(l + 2e)r 



(^-c) 2 + (^)-c) 2 |-ire( 



(5) 



The first two terms show that the inclusion locally promote an isotropic curvature c, and the third term distinguishes 
conformations for which the deviation of iJ-j; and iJ^ from c have identical or opposite signs. Note that it has been 
implicitely assumed that T > 0. 

In order to calculate the interaction energy between the two inclusions, we must minimize F m + U\ + U2 with 
respect to the shape of the membrane. We do this in two steps: (i) first, we minimize F m with respect to the 
membrane shape while setting prescribed values of H^, which yields F m - U1 (H^ , H a p)i (ii) next, we minimize F tot — 

(2) 



F, 



U2 with respect to H^I and H^) 



a/3 J -T ~U\ T w 2 wxuxi icopcv.j_.Aj M-L a0 aim j.± ^ 

The first step was performed in Ref. |6|. The result is 



(ttW tt(' 



1 



ft K* M -1 K. 



(G) 



In this expression, the subscript t means transpose, 



K t_/rr(l) rr(l) rrU) rr{2) rr(2) ff (2)s 



(7) 



is the 6x6 matrix: 



M 



/ Giin(O) 
G m2 (0) 
G 1122 (0) 
G mi (-r) 
G m2 (-r) 
\G U22 (-r) 



G m2 (0) 
G U22 (0) 

Gl222(0) 

G m2 (-r) 

Gn22(~ r ) 
Gi222( — r ) 



G U22 (Q) Gnu(r) G m2 (r) G n22 (r)\ 

Gi222(0) Gm 2 (r) Gn 22 (r) G1222O") 

G2222(0) Gii22(f) Gi222(r) G 2 222(r) 

G 1122 (-r) G im (0) G m2 (0) G 1122 (0) 

Gi222(~ r ) Gni2(0) Gn22(0) Gi222(0) 

G222 2 ( — r ) Gn22(0) Gi222(0) GtITI^)) 



(8) 
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in which r = rj r 2 , and G Q ^ 7< 5(r) = dadpd^dgG^r), where 

f d 2 q e iq r 

G(r) = J (2^p ^ 4 + C 2 g 2 + g 4 (l + ^? 2 + &89 4 ) (9) 

is the Green function associated with the elastic Hamiltonian. Note that no cutoff is required for the convergence of 
G a /3 7 s(0) thanks to terms of order q e and q 8 . The question of the high-wavevector cutoff will be raised later. 
To perform the second step of the minimization, we rewrite the interaction term as 

Ui + U 2 = ~ (K - C) 1 N (K - C) , (10) 

where 

r t_/ r ,(i) r (i) r (i) r (2) r (2) r (2)s ,,,, 
contains the curvatures attempted to be set by the inclusions, and N is the 6x6 block-diagonal matrix: 




N = 2e ® , (12) 



where £g> stands for the matrix direct product. 

The resulting total energy F to t = F nlln + U\ + U2 is quadratic in K, and its minimization yields 

Ftot.min = - y C l N [k M- 1 + T N] _1 NC + ~ C'NC. (13) 

After some algebra, this expression can be recast in the simpler form, similar to that of Eq. 10: 

Ftot,mh = ^C t (M + ^N- 1 )" 1 C. (14) 

Equation (|14|) is very general and corresponds to our central result: it describes the interaction between two pointwise 
inclusions to all orders in their inverse separation 1/r (as long as the Hamiltonian F m is fit to describe the membrane 
distortion energy at the scale r), for any value of the rigidities T and Te, and actually for isotropic or anisotropic 
inclusions, since the above derivation does not assume any form for the components of C. For anisotropic inclusions, 
C should contain the components of the curvature tensor favored by the inclusions . 

To proceed further, we must first determine the components of the matrix M. In the large r limit, i.e., r 3> 6g,6s 
(recall that b§ and bs are nanoscopic lengths), and assuming also r <C (i.e., taking £~ — > 00 and setting strictly 

l~ x = for the sake of simplicity), Eq. @ yields 

G(r)^C 1 + -i-r 2 ln(^). (15) 

167T G2 

The constants C\ and G2 depend on £ 7 , but they actually disappear in the calculation of G Q ^ 7 a(r), leaving l/(16ir)(x 2 + 
y 2 ) \n(x 2 + y 2 ) as the only relevant contribution in G(r). To calculate the components of the matrix M in the large 
r limit defined above, we set x = r and y = 0, yielding, Gnu = G1122 — — l/(47rr 2 ), G2222 — 3/(47rr 2 ), and 
G1112 = G1222 = 0. Next, we have 

/dq q 5 
16^ ^TWTToWTW) 55 w (16) 

where n a pjs contains the angular part of the integration on q. Simple calculations yield nun = "2222 = 3, rc.1122 = 1, 
and nni2 = 711222 = 0. For be rs 6s -C £j, ( a , as previously specified, the constant b is a nanometric length comparable 
to be or bs- Therefore, without loss of generality, we have (still in the large r limit): 



'3 1 

1 

1 3 




M - £2 I 11 1 I In ■,') Tf-r i <> 1 "I (',' u) CL") 
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To calculate the interaction Fi nt (r) between two identical isotropic inclusions, we set 

C t = (c,0,c,c,0,c), (18) 

and, using Eq. ifTIjl. we obtain: 

Fi n t(r) = ^ o-a+O -« • (19) 

W 2n 2 (k + 2Tb- 2 ) [ K + 4Tb- 2 (l + 2e)} 2 ^ \r e J y ' 

We now consider two limits: hard inclusions, corresponding to T/k ^> b 2 , and soft inclusions, corresponding to 
T/k <C b 2 (note that we keep e of order unity for the sake of simplicity). In the former case, the inclusions strongly 
set their preferred curvature; in the latter case, they only weakly tend to establish it. 

For hard inclusions, we let T go to infinity, and to leading order in 1/r, we obtain 

. , nb 2 c 2 fb\ A , , 

Fint.hardW ~ — — - . (20) 



64-7T 2 \r 

We thus recover the form obtained in Refs. 0, 01 , with n appearing in the numerator. Here the microscopic 
length b arising from the membrane Hamiltonian plays the role of the cutoff introduced in Refs. 0. [J Q . or the 
role of the radius of the disk modeling the inclusion in the case of Ref. 2} . 

For soft inclusions, we obtain to leading order in 1/r and at lowest-order in T (like in Ref. pj): 

F 3 c 2 2A 2 D 
* ^2^4 (1 + 2e ) = -^2^1 (1 + 2e )> ( 21 ) 

where k appears now to the second power in the denominator. This expression coincides exactly with the Eq. (37) of 
Ref. 0, in which e = was implicitely assumed. Contrary to the case of hard inclusions, the microscopic length b 
does not appear in the interaction, it would however reappear in the next order term C(r 4 ). We recall that 1 + 2e is 
positive, hence the interaction is always repulsive in the present case of isotropic inclusions. 

We have therefore shown that there is no fundamental discrepancy between the recent result of Marchenko and 
Misbah pj and the earlier results published in Refs. 0- 0- El S 01 ■ The result of Ref. El arrplies to soft inclusions, 
i.e., to inclusions that weakly tend to curve the membrane, while the result of Refs. |2j, y, U, [Ej, |6j is valid for hard 
inclusions that strongly enforce their preferred membrane curvature. 



III. THE ROLE OF THE CUT-OFF 



Let us now comment on the role of the high-wavevector cutoff. Of course, any well-defined clastic theory does involve 
a coarse-graining cutoff a -1 , and the elastic constants k(o), k§{cl), and na(a) depend on it. Thanks to the high-order 
bending elastic constants Ke and k$, there was no need to introduce a cutoff to calculate the interaction between the 
quadratic pointwise inclusions considered here. This is however somewhat arbitrary, since the microscopic length b 
arising from this procedure is expected to compare with the natural cutoff a. This simply means that reintrocuding 
the cutoff would merely redefine b, leaving its order of magnitude unchanged. It is therefore unimportant to specify 
whether the cutoff was introduced or not in the procedure defining the microscopic length b that appears in the 
interaction. 



IV. COARSE-GRAINING OF THE MEMBRANE-INCLUSION SYSTEM 



Describing an inclusion by a pointwise potential is of course a mathematical approximation hiding some physical mi- 
croscopic length-scale. One may actually choose the "resolution" of the pointwise description, i.e., the high-wavevector 
cutoff, by using a renormalization procedure. For instance, one may choose a "microscopic" cutoff comparable with 
the inverse inclusion size, or a much smaller cutoff corresponding to an elementary point size much larger than the 
inclusion. The question then naturally arises of how the constants of the pointwise quadratic potential will depend 
on the coarse-graining scale. 

Let us therefore consider a membrane described with an upper wavevector cutoff A = a -1 . The membrane hosts 
a single pointwise inclusion located at r = ri. We shall perform a coarse-graining of this system in order that the 
new wavevector cutoff becomes A' = (/ia) _1 < A. We assume the whole system to be in thermal equilibrium at the 
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temperature T, however, since all energies are quadratic, the renormalization will be athermal and we can set k^T = 1 
without loss of generality. The partition function can then be written as 

©ft exp(-F m [/i] -CMfc]), (22) 



where the inclusion's pointwise potential can be written as U x = |r(Ki - C 1 ) t N 1 (K 1 - Q), as in Eq. (|10l) . Here, Ni 

is the 3x3 uper-left block of the matrix N of Eq. 112fl . and = (H[\\ h[\ , i?22^)> Q = (c, 0, c). It is convenient to 
write exp(— U\) as the result of a Gaussian integral (Hubbard-Stratonovich transformation |{|): 

e X p(-Ui) oc Jd 3 $ exp -l&N^ 1 * + i& (c!- J -^K^e*'^ , (23) 

where Ki*(q) = f— q 2 h(q), — qiq2h(q), — <$h{cm. Then, going to the variables h(q), Eq. lf2*2^) becomes (up to a 
multiplicative constant): 

/ 1 \ r r r ^2„ „ 1 

(24) 



Z = JdHexpf-^&N^ + i&d^ Jj[dh(q)exp - f <? 4 |Mq)| 2 + ^%(q)e lqri 

Since the argument of the second exponential in the above equation is a quadratic function of /i(q), it is straightforward 
to integrate over the /i(q) with wavevectors in the range {^la)^ 1 < |q| < a" 1 . Once this integration is performed, the 
resulting integral over <!> is Gaussian, and yields, after some elementary algebra, 



Z = Jvh< cxp (-F m [h<] - Uf[h<}) , 



where h<(r) = (2tt)- 2 f q<1/{fia) d 2 q h(q)e^ r , and 



cff 



pcff 



(dad, 



C 



cff^eff 



r ctt e 



-9 a a Q /i < (ri)5 (3 9 / 3/i < (ri) 



■pert in 



nbl + 2T{\- M - 2 ) 



off 



e«bg + T (1 + 2e) (m 2 ~ 1) 
nblii 2 -4r(l + 2e) (m 2 -1) : 



(25) 

(26) 
(27) 

(28) 



where we have defined &o = 4v27ra. This set of exact equations is our second main result. It describes the effective 
quadratic coupling between the inclusion and the membrane in the coarse-grained pointwise description corresponding 
to the new cutoff (/^a) _1 . Note that the inclusion's preferred curvature Ci is unaffected by the coarse-graining 
procedure. This result may seem surprising, however one should realize that the actual curvature set by the inclusion 
will naturally be smaller in the coarse-grained description due to the fact that r eff < V. 

In principle, the interaction between two inclusions should be independent of the level of coarse-graining chosen to 
describe the inclusions (as long as their separation is larger than the length-scale of the cutoff) . One can indeed check 
that the interaction Eq. (jT§|l is unchanged under the substitutions b = b — > fxbo, T — > r cff and e — > e eff . 

Now, since the rigidities T and e are scale-dependent, an inclusion may be hard or soft depending on the coarse- 
graining level. Indeed, according to the invariance of the interaction mentioned above, the criterion for the hardness or 
softness of the inclusion is to compare T cS / k to (fibg) 2 . For fj, ^> 1, the "rigidity" T/(k&q), which becomes T eS / (n^?b^) 
after coarse-graining, is reduced by a factor ~/i 2 [l + 2F/(k6q)] 3> F Hence all inclusions effectively become soft at 
large coarse-graining scales. 

Is this in contradiction with the fact that hard inclusions have an interaction proportional to k while soft inclusions 
have an interaction inversely proportional to n 2 7 The answer is no, because the effective constants D and A introduced 
by Misbah et al. Q] actually depend on k. (This is another way to see that there is no discrepancy between the results 
of Ref. |l| and those of the previous works.) Indeed, starting from a hard inclusion (r S> k&o): we obtain, in the limit 
/i ^> 1 where all inclusions become soft, r cff ~ and e off ~ — j, while the preferred curvature Co is unchanged. 

With the correspondance D cS — —T cS and A cS — cor off , we obtain (A eS ) 2 D eS tx k 3 , hence Eq. (|21|l indeed yields an 
interaction proportional to k. 

We therefore conclude that the correct procedure to determine whether a given inclusions is fundamentally hard or 
soft is to estimate T at the coarse-graining level where the cutoff a compares with the inclusion size. The corresponding 
"microscopic" value T m must then, according to the discussion at the end of Sec. [HI be compared with Ka 2 (since b 
and a are expected to be of the same order of magnitude) . An inclusion will be fundamentally hard if T m 3> na 2 or 
fundamentally soft if T m <g; Ka 2 . 
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V. PROTEINS AS "HARD" INCLUSION 



We now argue that integral proteins of conical shape should be fundamentally considered as hard inclusions. First, 
it is very likely that the membrane elasticity is too weak to alter the conical shape of the protein, because of the 
strong dipolar interactions within the bundles of a-helices in the transmembrane domain. Beside, a significant change 
in the protein shape would most probably alter its function. We can therefore assume that the protein behaves as a 
rigid conical body and that the weakness of the curvature potential embodied in the parameter T m arises from the tilt 
degree of freedom of the lipids. Let us consider a conical protein having a circular section of radius R in the membrane 
plane and an aperture angle 28. An estimate of T m can be obtained by calculating the energy stored in the lipid tilt 
when one assumes that the membrane midplane remains perfectly flat. The latter is W rs BO 2 x 2ttR x £, where B 
is the elastic constant associated with the tilt and £ w (k/B) 1 / 2 the tilt relaxation length. With Co = 0/R, which 
corresponds to the membrane curvature that the protein attempts to set, this energy can be rewritten as W rs T>T m c^ 
with r m w 47Tfti? 3 /£. Hence we obtain 



r t? 3 

~ 4 
na 2 a 2 £ 



2 - 4 -727- ( 29 ) 



Interestingly, this quantity is independent of k and it grows as R 3 . With typically R ~ 6nm, a ~ £ ~ 3nm, we obtain 
T m /(Ka 2 ) « 100, which justifies treating conical proteins as a hard inclusion. 



VI. CONCLUSION 



We have calculated the induced interaction acting between pointwise membrane inclusions that locally interact 
up to quadratic order with the membrane curvature tensor. According to the rigidity of the potential constraining 
the local membrane curvature, inclusions can be viewed as hard or soft. The interaction between hard inclusions is 
proportional to k while the interaction between soft inclusions is inversely proportional to k 2 (both are also inversely 
proportional to the fourth power of the separation in the case of isotropic inclusions). This reconciliates the recent 
view of Marchenko and Misbah [ij with those of earlier works. 

To discriminate whether a given inclusion is fundamentally hard or soft, one should however consider the coupling 
constants defined at the microscopic level, i.e., defined in a pointwise description in which the elementary point-size 
(the cutoff of the model) compares with the inclusion size. Indeed, we have shown that coarse-graining the pointwise 
description to length-scales much larger that the inclusion's size transforms any hard inclusion into a soft one, however 
with coupling constants depending on k in such a way that the global interaction remains proportional to k if the 
inclusions was fundamentally hard. In the light of this analysis we have shown that transmembrane proteins should 
be considered as fundamentally hard inclusions. 
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